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Abstract. Recent advances in the study of electrodisintegration of ''He are presented and discussed. The 
pair-correlated hyperspherical harmonics method is used to calculate the initial and final state wave func- 
tions, with a realistic Hamiltonian consisting of the Argonne vis two-nucleon and Urbana IX three-nucleon 
interactions. The model for the nuclear current and charge operators retains one- and many-body contribu- 
tions. Particular attention is made in the construction of the two-body current operators arising from the 
momentum-dependent part of the two-nucleon interaction. Three-body current operators are also included 
so that the full current operator is strictly conserved. The present model for the nuclear current operator 
is tested comparing theoretical predictions and experimental data of pd radiative capture cross section and 
spin observables. 

PACS. 21.45.-l-v Few-body systems - 25.10.-l-s Nuclear reactions involving few-nucleon systems - 25.30.Dh 
Inelastic electron scattering to specific states - 25.30.Fj Inelastic electron scattering to continuum - 
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1 Introduction 

The theoretical study of the electromagnetic structure of 
few-body nuclei requires the knowledge of the nuclear wave 
functions and electromagnetic transition operators. In the 
case of processes involving two and/or three nucleons, it 
is possible to obtain very accurate bound- and scattering- 
state wave functions from realistic Hamiltonian models. 
Therefore, the different models for the nuclear electromag- 
netic current operator can be tested with the large variety 
of electromagnetic observables involving A=2 and 3 nu- 
clei. In particular, we concentrate our attention on the 
pd radiative capture and on the electrodisintegration of 
^He below and above deuteron breakup threshold (DBT). 
These processes have been extensively studied by sev- 
eral research groups (see Ref. ^ for a review). Most re- 
cently, the pd radiative capture and the ^He(e,e') reac- 
tions below DBT have been investigated by our group 
in Ref. I^j. The pair-correlated hyperspherical harmonics 
(PHH) method |3 has been used to calculate the A=3 
bound- and scattering-state wave functions from a real- 
istic Hamiltonian model consisting of the Argonne vis 
two-nucleon and Urbana IX three-nucleon -S interac- 
tions (AV18/UIX). The nuclear electromagnetic current 
operator included, in addition to the one-body convec- 
tion and spin-magnetization terms, also two-body contri- 
butions. These two-body terms were constructed following 



the method of Ref. 6 , with the goal of satisfying the cur- 
rent conservation relation (CCR) with the AV18. How- 
ever, within this method, only the dominant two-body 
terms, constructed from the momentum-independent part 
of the AV18, satisfy the CCR with this part of the in- 
teraction The two-body terms originated from the 
momentum-dependent part of the AVIS are not strictly 
conserved. In Ref. further transverse contribution, as- 
sociated with the p7r7 and lottj transition mechanisms and 
with the excitation of intermediate A resonances, were in- 
cluded. 

The main conclusions of Ref. j5j can be summarized 
as follows: the two-body contributions to the electromag- 
netic charge and current operators are crucial to achieve 
an overall good agreement between theory and experiment 
for all the observables under investigation. However, some 
discrepancies between theory and experiment have been 
found, in particular for the deuteron tensor polarization 
observables T20 and T21 of pd radiative capture at center- 
of-mass energy Ec.m. — 2 MeV, and, in Ref. also at 
Ec.m. — 3.33 MeV. In the analysis of Ref. [2] it has been 
suggested that these discrepancies might be due to the fact 
that the electromagnetic current operator satisfies only 
approximately the CCR with the nuclear Hamiltonian. In 
fact, when the T20 and T21 observables are calculated in 
the long- wavelength approximation (LWA), applying the 
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Siegert theorem, a quite good agreement with the data is 
obtained. 

In this work we present a new model for the nuclear 
current operator which exactly satisfies the CCR with the 
AV18/UIX Hamiltonian model. The model is then tested 
in the study of the pd radiative capture at the two pre- 
viously considered values of Ec.m., 2 and 3.33 MeV, and 
in the study of inclusive and exclusive electrodisintegra- 
tion of ■^He, both below and above DBT. The model for 
the nuclear electromagnetic current and charge operators 
is summarized in the following section. A detailed review 
will be given elsewhere jH)- Finally, the theoretical results 
are compared with the experimental data in Sec.|2| 



2 The nuclear electromagnetic charge and 
current operators 

The nuclear electromagnetic charge p(q) and current j(q) 
operators can be written as sums of one- and many-body 
terms that operate on the nucleon degrees of freedom. The 
one-body operators Pi{q) and ji(q) for particle i are de- 
rived from the non-relativistic reduction of the covariant 
single-nucleon current, by expanding in powers of 1/m, m 
being the nucleon mass ^l^ . The model commonly used |10| 
for the two-body charge operators includes the tt-, p-, and 
w-meson exchange terms with both isoscalar and isovector 
components, as well as the (isoscalar) pnj and (isovector) 
charge transition couplings. At moderate values of 
momentum transfer (g < 5 fm"-'^), the 7r-meson exchange 
charge operator has been found to give the dominant two- 
body contribution 

The electromagnetic current operator must satisfy the 
CCR, written as 



q.j(q) = [i7,p(q)] 



(1) 



where the nuclear Hamiltonian H is taken to consist of 
two- and three-body interactions, denoted as Vij and Vijk 
respectively. To lowest order in 1/m, Eq. separates into 



[2m 

q- j-y(q) = b*j,Pi(q) + Pj(q)] , 



(2) 
(3) 



and similarly for the three-body current jyfe(q). We have 
neglected the two-body terms in p(q), which are of or- 
der Ijiv?. The one-body current is easily shown to satisfy 
Eq. |(2J). To construct the two-body current, it is useful 
to adopt the classification scheme of Ref. , and sepa- 
rate the current ]ij (q) into model-independent (MI) and 
model-dependent (MD) parts. The MI two-body current 
has a longitudinal component and is constructed so as to 
satisfy the CCR of Eq. (jSJ, while the MD two-body cur- 
rent is purely transverse and therefore is un-constrained 
by the CCR. The latter is taken to consist of the isoscalar 
/37r7 and isovector wk^ transition currents, as well as the 
isovector current associated with excitation of intermedi- 
ate A resonances as in Ref. 



The MI two-body currents arising from the momentum- 
independent terms of the AV18 two-nucleon interaction 
have been constructed following the standard procedure of 
Ref. [HI , which will be hereafter quoted as meson-exchange 
(ME) scheme. It can be shown that these two-body current 
operators satisfy exactly the CCR with the first six opera- 
tors of the AV18. The two-body currents arising from the 
spin-orbit components of the AVIS could be constructed 
using again ME mechanisms |13| . but the resulting cur- 
rents turn out to be not strictly conserved. The same 
can be said of those currents deriving from the quadratic 
momentum-dependent components of the AVIS, if obtained, 
as in Ref. j2j, by gauging only the momentum operators, 
but ignoring the implicit momentum dependence which 
comes through the isospin exchange operator (see below) . 
Since our goal is to construct MI two-body currents which 
satisfy exactly the CCR of Eq. the currents arising 
from the momentum-dependent terms of the AVIS in- 
teraction have been obtained following the procedure of 
Ref. 14). which will be quoted as minimal-substitution 
(MS) scheme. The main idea of this procedure, reviewed 
and extended in Ref. |H] , is that the isospin operator • Tj , 
is formally equivalent to an implicit momentum depen- 
dence jl4) . In fact, Ti ■ Tj can be expressed in terms of the 
space-exchange operator {Pij) using the formula 



-1-(1 + . 



(4) 



valid when operating on antisymmetric wave functions. 
The operator Pij is defined as Pij = e''j''^'+'''J'^3 , where 
the V-operators do not act on the vectors r.^ = — Vj = 
—Vji in the exponential. In the presence of an electro- 
magnetic field, minimal substitution is performed both in 
the momentum dependent terms of the two-nucleon inter- 
action and in the space-exchange operator Pij of Eq. I^J , 
and the resulting current operators are then obtained with 
standard procedures |MT^ . Explicit formulas can be found 
in Ref. Here we only quote the result for the currents 
arising from the momentum-independent part of the AVIS 
interaction (u° ): 



7,. 



(5) 

where ti is the nucleon charge operator, and ds [ds') is the 
infinitesimal step on the generic path ^ij (7^^) that goes 
from position i (j) to position j (i). Two observations are 
in order: (i) with a particular choice of the integration path 
it is possible to re-obtain the two-body current operators 
calculated, within the ME scheme, from the momentum- 
independent part of the AVIS 9 ; (ii) in the limit q — » 
0, the current operator jij(q) of Eq. lO becomes path- 
independent and unique. To simplify the calculation, the 
integration paths 7^ and 7^^ have been chosen to be linear. 

Both the ME and the MS schemes can be general- 
ized to calculate the three-body current operators induced 
by the three- nucleon interaction (TNI) Vijk- Here, these 
three-body currents have been constructed within the ME 
scheme to satisfy the CCR with the Urbana-IX TNI 0. 
However, the same procedure can be applied to other 
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Fig. 1. Deuteron tensor polarization observables T20 and T21 
for pd radiative capture at Ec.m.~ 2 MeV, obtained with the 
AVI 8 Hamiltonian models. See text for explanations of the 
different curves. The experimental data are from Ref. |19|. 

models of TNI's, such as the Tucson-Melbourne |E1 and 
Brazil ^H] niodels. Details of the calculation can be found 
in Ref. 0. 

In summary, the present model for the many-body cur- 
rent operator retains the two-body terms obtained within 
the ME scheme from the momentum-independent part of 
the AV18, those ones obtained within the MS scheme 
from the momentum-dependent part of the AV18, the 
MD terms quoted above, and the three-body terms ob- 
tained within the ME scheme from the UIX. Thus, the 
full current operator satisfies exactly the CCR with the 
AV18/UIX nuclear Hamiltonian. In contrast, the model 
of Ref. retains only two-body currents, all of them ob- 
tained within the ME scheme. 




Fig. 2. Differential cross section, proton vector analyzing 
power, and the four deuteron tensor polarization observables 
for pd radiative capture at Ec.m. =3.33 MeV, obtained with 
the AV18/UIX Hamiltonian model. See text for the explana- 
tions of the different curves. The experimental data are from 
Ref. 1201. 



3 Results 

In the present section we compare the theoretical predic- 
tion and experimental data for two sets of "^He electrodis- 
integration observables: (i) the longitudinal and transverse 
response functions Rl and Rt for ■^IIe(e, e') at momentum 
transfer values (7=0.88, 1.64 and 2.47 fm~^ and excitation 
energies (Ex) from two-body threshold up to 20 MeV; (ii) 
the differential cross section of the '^IIe(e, e'd)p reaction 
as function of the missing momentum in (g, w)-constant 
kinematics, at beam energies of 370 and 576 MeV and q 
values of 412, 504 and 604 MeV/c. Data are respectively 
from Refs. |I7| and [TH|. 

However, prehminarly, we report some results for the 
pd radiative capture reaction, to show that the new model 
for the many-body nuclear current operators resolves some 
of the discrepancies between theory and experiment of 
Ref. m. The T20 and T21 observables at Ec.m.='2 MeV us- 
ing the AV18 two-nucleon interaction are shown in Fig.^ 
The dotted curves are obtained including only the one- 
body current contributions, the dashed curves are ob- 
tained with the model for the nuclear current operator 
of Ref. 121, and the dotted-dashed curves are obtained in 
the long- wavelength-approximation (LWA), applying the 



Siegert theorem. Finally, the solid curves are obtained in- 
cluding the one- and two-body contributions described in 
Sec.|5J necessary to satisfy the CCR with the AV18 nuclear 
Hamiltonian. Data are from Ref. jEl- By inspection of 
Fig. n] we observe that there is a good agreement between 
the experimental data, the LWA and the present "full" 
results, while the results of Ref. are in disagreement, 
as expected, with both the LWA and experimental results. 
The excellent agreement between the LWA and the "full" 
results is a consequence of the fact that the "full" nuclear 
electromagnetic current operator satisfies the CCR with 
the AV18 nuclear Hamiltonian. 

As an example of the degree of agreement which has 
been reached between the present calculation and experi- 
mental data for the pd radiative capture, we show in Fig.|2 
the theoretical predictions obtained with the AV18/UIX 
Hamiltonian model at £^c.m. =3.33 MeV. Data are from 
Ref. 1201 • The dashed, dotted-dashed and solid curves cor- 
respond to the calculation with one-body only, with one- 
and two-body, and with one-, two- and three-body cur- 
rents. An overall nice description has been reached for 
all the observables, with the only exception of the iTn 
deuteron polarization observable at small angles. Also, 
some small three-body currents effects are noticeable, es- 
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Fig. 3. Longitudinal and transverse response functions of '^He, 
obtained with the AV18/UIX Hamiltonian model. The exper- 
imental data are from Ref. 1171 . The vertical line represents 
the ppn breakup threshold. See text for explanations of the 
different curves. 



pecially in the T20 and T21 deuteron tensor observables, 
which is an indication of the fact that if a Hamiltonian 
model with two- and three-nucleon interactions is used, 
then the model for the nuclear current operator should 
include the corresponding two- and three-body contribu- 
tions. 

The theoretical predictions for the longitudinal and 
transverse response functions Rl and Rt for g=0.88, 1.64 
and 2.47 fm^i and 5 MeV < < 20 MeV are shown 
in Fig. 131 The vertical lines indicate the DBT. Although 
the calculation is extended above this threshold, explicit 
three-body breakup channel contributions have not yet 
been included. The dashed, dotted-dashed and solid lines 
are obtained with only one-body current and charge op- 
erators, with the one- and two-body operators of Ref. 2 , 
and with the one- and many-body transition operators 
presented in Sec. |2| Note that in the case of the charge 
operator, the model of Ref. [2] and the present one coin- 
cide. From inspection of the figure, we can observe that 
there is no significant difference between the results ob- 
tained using the model for the nuclear current operator of 
Ref. PI and the present one. To draw any conclusion in the 
comparison between theory and experiment, a complete 
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Fig. 4. Differential cross section of the ^He(e, e'd)p reaction 
obtained with the AV18/UIX Hamiltonian model as function 
of the missing momentum at three q values and two beam en- 
ergies. The data are from Ref. |18j . See text for the explanation 
of the different curves. 



calculation which includes also the ppn channel should be 
performed. 

The theoretical predictions for the differential cross 
section of the ''He(e, e'd)p reaction as function of the miss- 
ing momentum at beam energies of 370 and 576 MeV and 
q values of 412, 504 and 604 MeV/c are compared with 
the data of Ref. TH* in Fig. El The dotted-dashed lines 
correspond to the plane-wave impulse approximation re- 
sults. When final-state-interaction effects are included and 
the pd final wave function is calculated with the PHH 
technique using the AV18/UIX Hamiltonian model, the 
dashed and the solid lines are obtained, depending if the 
one-body only or the one- and many-body contributions 
to the nuclear transition operators are retained. Here, no 
comparison is shown between the old model of Ref. [2] 
and the present one, since no significant differences have 
been seen. Note that this is the first calculation for the 
'^He(e, e'd)p reaction above DBT which uses the PHH tech- 
nique to calculate the bound- and scattering-state wave 
functions including TNI and Coulomb force effects in the 
final-state interaction. By inspection of Fig. ^ we can 
conclude that there is an overall nice agreement between 
theory and experiment, although the data at low missing 
momentum are overestimated by the theory. This was al- 
ready observed in Ref. 18 , where the data were compared 
with a Faddeev calculation, with no Coulomb and three- 
nucleon interaction. However, we have verified that the 
inclusion of the UIX three-nucleon interaction improves 
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Fig. 5. Differential cross section of the ^He(e, e'd)p reaction 
as function of the missing momentum at q = 412 MeV/c 
and beam energy of 370 MeV and 576 MeV. The dashed and 
solid lines correspond to the calculation performed with the 
AV18 and with the AV18/UIX Ifamiltonian models, respec- 
tively. Data are from Ref. |18|. 



the description of the data at zero missing momentum. 
This can be seen in Fig. |S1 where the differential cross 
section at q = 412 MeV/c and beam energy of 370 and 
576 MeV is calculated vifith the AV18 (dashed curves) and 
the AV18/UIX (solid curves) Hamiltonian models. Only 
in this second case the model for the nuclear current op- 
erator includes three-body contributions, so that for each 
given Hamiltonian model the CCR is satisfied. 



4 Summary and Outlook 

We have reported on new calculations of ^He(e, e') longi- 
tudinal and transverse response functions for three values 
of the momentum transfer and excitation energies from 
two-body threshold up to 20 MeV, and for '^He(e, e'd)p 
differential cross section as function of the missing mo- 
mentum in {q, ti;)-constant kinematics, at two beam ener- 
gies and three q values. These calculations use accurate 
bound and scattering state wave functions obtained with 
the PHH method from the Argonne uig two-nucleon and 
Urbana IX three-nucleon interactions. The model for the 
electromagnetic charge operator includes one- and two- 
body components, while the model for the electromag- 
netic current operator includes one-, two- and three-body 
components, constructed so as to satisfy exactly the CCR 
with the given Hamiltonian model. The model for the nu- 
clear current operator has been tested calculating the pd 
radiative capture cross section and spin observables. In 
particular, we have shown that the experimental T20 and 
T21 deuteron tensor observables are nicely reproduced by 
the present calculation. A systematic comparison between 
theory and experiment for the pd radiative capture in a 
wide range of Ec.m. is currently underway |^. 

The '^He electrodisintegration observables considered 
in the present work are very sensitive to the many-body 
contributions of the nuclear transition operators. However, 
no significant differences have been observed between the 
calculation performed with the present model of the nu- 
clear current operator and the one of Ref. 'S' . In the case of 



the '^He(e, e') reaction, theoretical predictions and experi- 
mental data cannot be directly compared above deuteron 
breakup threshold, since the (energetically open) three- 
body breakup channel contributions have yet to be cal- 
culated explicitly. In the case of the ■^He(e, e'd)p reaction, 
instead, the experimental differential cross section is fairly 
reproduced in the whole range of missing momentum for 
all different kinematics, when final-state interactions ef- 
fects are considered. Further investigations for the •^He 
electrodisintegration in a wider energy range are vigor- 
ously being pursued. 

The work of R.S. was supported by the U.S. DOE Con- 
tract No. DE-AC05-84ER40150, under which the South- 
eastern Universities Research Association (SURA) oper- 
ates the Thomas Jefferson National Accelarator Facility. 
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1 Introduction 

The theoretical study of the electromagnetic structure of 
few-body nuclei requires the knowledge of the nuclear wave 
functions and electromagnetic transition operators. In the 
case of processes involving two and/or three nucleons, it 
is possible to obtain very accurate bound- and scattering- 
state wave functions from reahstic Hamiltonian models. 
Therefore, the different models for the nuclear electromag- 
netic current operator can be tested with the large variety 
of electromagnetic observables involving A=2 and 3 nu- 
clei. In particular, we concentrate our attention on the 
pd radiative capture and on the electrodisintegration of 
^He below and above deuteron breakup threshold (DBT). 
These processes have been extensively studied by sev- 
eral research groups (see Ref. [1] for a review). Most re- 
cently, the pd radiative capture and the ^He(e,e') reac- 
tions below DBT have been investigated by our group 
in Ref. [2]. The pair-correlated hyperspherical harmonics 
(PHH) method [3] has been used to calculate the A=3 
bound- and scattering-state wave functions from a real- 
istic Hamiltonian model consisting of the Argonne wig 
two-nucleon [4] and Urbana IX three-nucleon [5] interac- 
tions (AV18/UIX). The nuclear electromagnetic current 
operator included, in addition to the one-body convec- 
tion and spin-magnetization terms, also two-body contri- 
butions. These two-body terms were constructed following 



the method of Ref. [6] , with the goal of satisfying the cur- 
rent conservation relation (CCR) with the AV18. How- 
ever, within this method, only the dominant two-body 
terms, constructed from the momentum- independent part 
of the AV18, satisfy the CCR with this part of the in- 
teraction [7]. The two-body terms originated from the 
momentum-dependent part of the AVIS are not strictly 
conserved. In Ref. [2], further transverse contribution, as- 
sociated with the p7r7 and cott^ transition mechanisms and 
with the excitation of intermediate A resonances, were in- 
cluded. 

The main conclusions of Ref. [2] can be summarized 
as follows: the two-body contributions to the electromag- 
netic charge and current operators are crucial to achieve 
an overall good agreement between theory and experiment 
for all the observables under investigation. However, some 
discrepancies between theory and experiment have been 
found, in particular for the deuteron tensor polarization 
observables T20 and T21 of pd radiative capture at center- 
of-mass energy Ecm. = 2 MeV, and, in Ref. [8], also at 
Ec.m. = 3.33 MeV. In the analysis of Ref. [2] it has been 
suggested that these discrepancies might be due to the fact 
that the electromagnetic current operator satisfies only 
approximately the CCR with the nuclear Hamiltonian. In 
fact, when the T20 and T21 observables are calculated in 
the long- wavelength approximation (LWA), applying the 
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Siegert theorem, a quite good agreement with the data is 
obtained. 

In this work we present a new model for the nuclear 
current operator which exactly satisfies the CCR with the 
AV18/UIX Hamiltonian model. The model is then tested 
in the study of the pd radiative capture at the two pre- 
viously considered values of Ec.m., 2 and 3.33 MeV, and 
in the study of inclusive and exclusive electrodisintegra- 
tion of ^He, both below and above DBT. The model for 
the nuclear electromagnetic current and charge operators 
is summarized in the following section. A detailed review 
will be given elsewhere [9]. Finally, the theoretical results 
are compared with the experimental data in Sec. 3. 



2 The nuclear electromagnetic charge and 
current operators 

The nuclear electromagnetic charge p{q) and current j(q) 
operators can be written as sums of one- and many-body 
terms that operate on the nucleon degrees of freedom. The 
one-body operators Pi{q) and j4(q) for particle i are de- 
rived from the non-relativistic reduction of the covariant 
single-nucleon current, by expanding in powers of 1/m, m 
being the nucleon mass [1]. The model commonly used [10] 
for the two-body charge operators includes the tt-, p-, and 
w-meson exchange terms with both isoscalar and isovector 
components, as well as the (isoscalar) pirj and (isovector) 
coTTj charge transition couplings. At moderate values of 
momentum transfer (g < 5 fm^^), the 7r-meson exchange 
charge operator has been found to give the dominant two- 
body contribution [11]. 

The electromagnetic current operator must satisfy the 
CCR, written as 



(1) 



where the nuclear Hamiltonian H is taken to consist of 
two- and three-body interactions, denoted as Vij and Vtjk 
respectively. To lowest order in 1/m, Eq. (1) separates into 



qji(q) 



2m 



.p»(q) 



q-jij(q) = [vij,Pii(i) +Pj(q)] , 



(2) 
(3) 



and similarly for the three-body current jyj; (q). We have 
neglected the two-body terms in p(q), which are of or- 
der 1/m^. The one-body current is easily shown to satisfy 
Eq. (2). To construct the two-body current, it is useful 
to adopt the classification scheme of Ref. [12], and sepa- 
rate the current jy (q) into model-independent (MI) and 
model-dependent (MD) parts. The MI two-body current 
has a longitudinal component and is constructed so as to 
satisfy the CCR of Eq. (3), while the MD two-body cur- 
rent is purely transverse and therefore is un-constraincd 
by the CCR. The latter is taken to consist of the isoscalar 
pTT'f and isovector a;7r7 transition currents, as well as the 
isovector current associated with excitation of intermedi- 
ate A resonances as in Ref. [2]. 



The MI two-body currents arising from the momentum- 
independent terms of the AV18 two-nucleon interaction 
have been constructed following the standard procedure of 
Ref. [6] , which will be hereafter quoted as meson-exchange 
(ME) scheme. It can be shown that these two-body current 
operators satisfy exactly the CCR with the first six opera- 
tors of the AV18. The two-body currents arising from the 
spin-orbit components of the AVIS could be constructed 
using again ME mechanisms [13], but the resulting cur- 
rents turn out to be not strictly conserved. The same 
can be said of those currents deriving from the quadratic 
momentum-dependent components of the AV18, if obtained, 
as in Ref. [2], by gauging only the momentum operators, 
but ignoring the implicit momentum dependence which 
comes through the isospin exchange operator (see below) . 
Since our goal is to construct MI two-body currents which 
satisfy exactly the CCR of Eq. (3), the currents arising 
from the momentum-dependent terms of the AV18 in- 
teraction have been obtained following the procedure of 
Ref. [14], which will be quoted as minimal-substitution 
(MS) scheme. The main idea of this procedure, reviewed 
and extended in Ref. [9] , is that the isospin operator r, ■ Tj , 
is formally equivalent to an implicit momentum depen- 
dence [14]. In fact, Ti ■ Tj can be expressed in terms of the 
space-exchange operator (Pjj) using the formula 



-1- (l + at ■ (Tj)Pij 



(4) 



valid when operating on antisymmetric wave functions. 
The operator is defined as Pij = e'"j''^'+'"'j '^j , where 
the V-operators do not act on the vectors Vij = Vi — Vj = 
—Tji in the exponential. In the presence of an electro- 
magnetic field, minimal substitution is performed both in 
the momentum dependent terms of the two-nucleon inter- 
action and in the space-exchange operator Py of Eq. (4) , 
and the resulting current operators are then obtained with 
standard procedures [9, 14]. Explicit formulas can be found 
in Ref. [9]. Here we only quote the result for the currents 
arising from the momentum-independent part of the AV18 
interaction (vfj): 

jy(q) = iv°Jei I dse'^-^ + Ej J ds' e'^A {l + n-Tj) , 

(5) 

where ej is the nucleon charge operator, and ds {ds') is the 
infinitesimal step on the generic path (7^^) that goes 
from position i (j) to position j (i). Two observations are 
in order: (i) with a particular choice of the integration path 
it is possible to re-obtain the two-body current operators 
calculated, within the ME scheme, from the momentum- 
independent part of the AV18 [9]; (ii) in the limit q 
0, the current operator jy(q) of Eq. (5) becomes path- 
independent and unique. To simplify the calculation, the 
integration paths jij and 7^^ have been chosen to be linear. 

Both the ME and the MS schemes can be general- 
ized to calculate the three-body current operators induced 
by the three-nuclcon interaction (TNI) Vijk- Here, these 
three-body currents have been constructed within the ME 
scheme to satisfy the CCR with the Urbana-IX TNI [5]. 
However, the same procedure can be applied to other 
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Fig. 1. Deuteron tensor polarization observables T20 and T21 
for pd radiative capture at Ec.m.= 2 MeV, obtained with the 
AV18 Hamiltonian models. See text for explanations of the 
different curves. The experimental data are from Ref. [19]. 

models of TNFs, such as the Tucson-Melbourne [15] and 
Brazil [16] models. Details of the calculation can be found 
in Ref. [9]. 

In summary, the present model for the many-body cur- 
rent operator retains the two-body terms obtained within 
the ME scheme from the momentum-independent part of 
the AV18, those ones obtained within the MS scheme 
from the momentum-dependent part of the AV18, the 
MD terms quoted above, and the three-body terms ob- 
tained within the ME schemo from the UIX. Thus, the 
full current operator satisfies exactly the CCR with the 
AV18/UIX nuclear Hamiltonian. In contrast, the model 
of Ref. [2] retains only two-body currents, all of them ob- 
tained within the ME scheme. 




8c Weg] [deg] 




30 60 90 120 150 180 ■ 30 60 90 120 150 180 
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Fig. 2. Differential cross section, proton vector analyzing 
power, and the four deuteron tensor polarization observables 
for pd radiative capture at Ec.m. =3.33 MeV, obtained with 
the AV18/UIX Hamiltonian model. See text for the explana- 
tions of the different curves. The experimental data are from 
Ref. [20]. 



3 Results 

In the present section we compare the theoretical predic- 
tion and experimental data for two sets of ^He electrodis- 
integration observables: (i) the longitudinal and transverse 
response functions and Rt for ^He(e, e') at momentum 
transfer values g=0.88, 1.64 and 2.47 fm~^ and excitation 
energies {E^) from two-body threshold up to 20 MeV; (ii) 
the differential cross section of the ^He(e,e'(i)p reaction 
as function of the missing momentum in (g, a;)-constant 
kinematics, at beam energies of 370 and 576 MeV and q 
values of 412, 504 and 604 MeV/c. Data are respectively 
from Refs. [17] and [18]. 

However, preliminarly, we report some results for the 
pd radiative capture reaction, to show that the new model 
for the many-body nuclear current operators resolves some 
of the discrepancies between theory and experiment of 
Ref. [2]. The T20 and T21 observables at £c.m.=2 MeV us- 
ing the AV18 two-nucleon interaction are shown in Fig. 1. 
The dotted curves are obtained including only the one- 
body current contributions, the dashed curves are ob- 
tained with the model for the nuclear current operator 
of Ref. [2], and the dotted-dashed curves are obtained in 
the long- wavelength-approximation (LWA), applying the 



Siegert theorem. Finahy, the sohd curves are obtained in- 
cluding the one- and two-body contributions described in 
Sec. 2, necessary to satisfy the CCR with the AV18 nuclear 
Hamiltonian. Data are from Ref. [19]. By inspection of 
Fig. 1, we observe that there is a good agreement between 
the experimental data, the LWA and the present "full" 
results, while the results of Ref. [2] are in disagreement, 
as expected, with both the LWA and experimental results. 
The excellent agreement between the LWA and the "full" 
results is a consequence of the fact that the "full" nuclear 
electromagnetic current operator satisfies the CCR with 
the AVI 8 nuclear Hamiltonian. 

As an example of the degree of agreement which has 
been reached between the present calculation and experi- 
mental data for the pd radiative capture, we show in Fig. 2 
the theoretical predictions obtained with the AV18/UIX 
Hamiltonian model at Sc „j,=3.33 MeV. Data are from 
Ref. [20]. The dashed, dotted-dashed and solid curves cor- 
respond to the calculation with one-body only, with one- 
and two-body, and with one-, two- and three-body cur- 
rents. An overall nice description has been reached for 
all the observables, with the only exception of the iTn 
deuteron polarization observable at small angles. Also, 
some small three-body currents effects are noticeable, es- 
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Fig. 3. Longitudinal and transverse response functions of ^He, 
obtained with the AV18/UIX Hamiltonian model. The exper- 
imental data are from Ref. [17]. The vertical line represents 
the ppn breakup threshold. See text for explanations of the 
different curves. 



pecially in the T20 and T21 deuteron tensor observables, 
which is an indication of the fact that if a Hamiltonian 
model with two- and three-nucleon interactions is used, 
then the model for the nuclear current operator should 
include the corresponding two- and three-body contribu- 
tions. 

The theoretical predictions for the longitudinal and 
transverse response functions and Rt for g=0.88, 1.64 
and 2.47 fm'^ and 5 MeV < < 20 MeV are shown 
in Fig. 3. The vertical lines indicate the DBT. Although 
the calculation is extended above this threshold, explicit 
three-body breakup channel contributions have not yet 
been included. The dashed, dotted-dashed and solid Hues 
are obtained with only one-body current and charge op- 
erators, with the one- and two-body operators of Ref. [2], 
and with the one- and many-body transition operators 
presented in Sec. 2. Note that in the case of the charge 
operator, the model of Ref. [2] and the present one coin- 
cide. From inspection of the figure, we can observe that 
there is no significant difference between the results ob- 
tained using the model for the nuclear current operator of 
Ref. [2] and the present one. To draw any conclusion in the 
comparison between theory and experiment, a complete 
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Fig. 4. Differential cross section of the ^He(e, e'd)p reaction 
obtained with the AV18/UIX Hamiltonian model as function 
of the missing momentum at three q values and two beam en- 
ergies. The data are from Ref. [18]. See text for the explanation 
of the different curves. 



calculation which includes also the ppn channel should be 
performed. 

The theoretical predictions for the differential cross 
section of the ^He(e, e'd}p reaction as function of the miss- 
ing momentum at beam energies of 370 and 576 MeV and 
q values of 412, 504 and 604 MeV/c are compared with 
the data of Ref. [18] in Fig. 4. The dotted-dashed lines 
correspond to the plane-wave impulse approximation re- 
sults. When final-state-interaction effects are included and 
the pd final wave function is calculated with the PHH 
technique using the AV18/UIX Hamiltonian model, the 
dashed and the sohd lines are obtained, depending if the 
one-body only or the one- and many-body contributions 
to the nuclear transition operators are retained. Here, no 
comparison is shown between the old model of Ref. [2] 
and the present one, since no significant differences have 
been seen. Note that this is the first calculation for the 
^He(e, e'd)p reaction above DBT which uses the PHH tech- 
nique to calculate the bound- and scattering-state wave 
functions including TNI and Coulomb force effects in the 
final-state interaction. By inspection of Fig. 4, we can 
conclude that there is an overall nice agreement between 
theory and experiment, although the data at low missing 
momentum are overestimated by the theory. This was al- 
ready observed in Ref. [18], where the data were compared 
with a Faddeev calculation, with no Coulomb and three- 
nucleon interaction. However, we have verified that the 
inclusion of the UIX three-nucleon interaction improves 
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Fig. 5. Differential cross section of the ^He(e, e'd)p reaction 
as function of the missing momentum at q = 412 MeV/c 
and beam energy of 370 MeV and 576 MeV. The dashed and 
solid lines correspond to the calculation performed with the 
AV18 and with the AV18/UIX Hamiltonian models, respec- 
tively. Data are from Ref. [18]. 

the description of the data at zero missing momentum. 
This can be seen in Fig. 5, where the diflerential cross 
section at g = 412 MeV/c and beam energy of 370 and 
576 MeV is calculated with the AV18 (dashed curves) and 
the AV18/UIX (soHd curves) Hamiltonian models. Only 
in this second case the model for the nuclear current op- 
erator includes three-body contributions, so that for each 
given Hamiltonian model the CCR is satisfied. 

4 Summary and Outlook 

We have reported on new calculations of ^He(e,e') longi- 
tudinal and transverse response functions for three values 
of the momentum transfer and excitation energies from 
two-body threshold up to 20 MeV, and for •^He(e, e'd)p 
differential cross section as function of the missing mo- 
mentum in {q, w)-constant kinematics, at two beam ener- 
gies and three q values. These calculations use accurate 
bound and scattering state wave functions obtained with 
the PHH method from the Argonne vis two-nucleon and 
Urbana IX three-nucleon interactions. The model for the 
electromagnetic charge operator includes one- and two- 
body components, while the model for the electromag- 
netic current operator includes one-, two- and three-body 
components, constructed so as to satisfy exactly the CCR 
with the given Hamiltonian model. The model for the nu- 
clear current operator has been tested calculating the pd 
radiative capture cross section and spin observables. In 
particular, we have shown that the experimental T20 and 
T21 deuteron tensor observables are nicely reproduced by 
the present calculation. A systematic comparison between 
theory and experiment for the pd radiative capture in a 
wide range of Ec.m. is currently underway [9]. 

The ^He electrodisintegration observables considered 
in the present work are very sensitive to the many-body 
contributions of the nuclear transition operators. However, 
no significant differences have been observed between the 
calculation performed with the present model of the nu- 
clear current operator and the one of Ref. [2] . In the case of 



the ^He(e, e') reaction, theoretical predictions and experi- 
mental data cannot be directly compared above deuteron 
breakup threshold, since the (energetically open) three- 
body breakup channel contributions have yet to be cal- 
culated explicitly. In the case of the ^He(e,e'd)p reaction, 
instead, the experimental differential cross section is fairly 
reproduced in the whole range of missing momentum for 
all different kinematics, when final-state interactions ef- 
fects are considered. Further investigations for the ^He 
electrodisintegration in a wider energy range are vigor- 
ously being pursued. 
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